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This paper continue earlier investigations on the decay of Burgers turbulence in one dimension 
from Gaussian random initial conditions of the power-law spectral type Eo{k) ~ |fc|". Depending 
on the power n, different characteristic regions are distinguished. The main focus of this paper is 
to delineate the regions in wave-number k and time t in which self-similarity can (and cannot) be 
observed, taking into account small-fc and large-fc cutoffs. The evolution of the spectrum can be 
inferred using physical arguments describing the competition between the initial spectrum and the 
new frequencies generated by the dynamics. For large wavenumbers, we always have k~^ region, 
associated to the shocks. When n is less than one, the large-scale part of the spectrum is preserved in 
time and the global evolution is self-similar, so that scaling arguments perfectly predict the behavior 
in time of the energy and of the integral scale. If n is larger than two, the spectrum tends for long 
times to a universal scaling form independent of the initial conditions, with universal behavior k^ 
at small wavenumbers. In the interval 2 < n the leading behaviour is self-similar, independent of 
n and with universal behavior k^ at small wavenumber. When 1 < n < 2, the spectrum has three 
scaling regions : first, a Ifcl" region at very small ks with a time-independent constant, second, a k^ 
region at intermediate wavenumbers, finally, the usual k~^ region. In the remaining interval, n < —3 
the small-fc cutoff dominates, and n also plays no role. We find also (numerically) the subleading 
term ~ k^ in the evolution of the spectrum in the interval —3 < n < 1. High-resolution numerical 
simulations have been performed confirming both scaling predictions and analytical asymptotic 
theory. 
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I. INTRODUCTION 



We study here Burgers equation 



dv 



dv 
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dx dx^ 



(1) 



in the limit of vanishing coefRcient v. First introduced 
by J.M. Burgers as a model of hydrodynamic turbulence, 
this equation arises in many situations in physics, see 
1^ m IE IS for classical (and more recent) mono- 
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graphs. It is fair to say that one of the main interest 
in Burgers equation over the last decade has been as a 
model for structure formation in the earlyuniverse within 
the so-called adhesion approximation 0)M0|- The Hopf- 
Cole transformation, to which we will return below, has 
been developped into a powerful tool to elucidate the sta- 
tistical properties of solutions to Burgers equations with 
random initial conditions of cosmological type [9l ll0lllll |. 
If a random force is added to the right-hand side of Q), 
the resulting KPZ equation is one of the most important 
models of e.g. surface growth [T3.ll3l[T^ . 

Investigations of Burgers turbulence have a long pre- 
history, started already by Burgers (1974), who was 
mainly concerned with white noise initial conditions. But 
nevertheless only recently jl^ the exact statistical prop- 
erties of the Burgers equation for the case n = were 
found. The case of fractal Brownian motion for the po- 
tential or for initial velocity is much more complicated 
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Burgers equation describes two principal effects in- 
herent in any turbulence |17| : the nonlinear redistribu- 
tion of energy over the spectrum and the action of vis- 
cosity in small-scale regions. Except for the direct phys- 
ical application, Burgers equation is hence also of great 
interest to test theories and models of fully developped 
turbulence. This paper follows that tradition. In an ear- 
lier contribution 18] we showed how self-similarity argu- 
ments, going back to Kolmogorov and Loitsyanski 
pof can be be disproved, in Burgers equation, for a class 
of initial conditions. A similar result was later arrived 
at by Eyink and Thompson for the Navier-Stokes equa- 
tion [21|, within an eddy-damped, quasi-normal Marko- 
vian (EDQNM) scheme. In this paper, we will discuss in 
greater detail how the self-similar (and not self-similar) 
regimes are realized with initial conditions that are only 
self-similar over a finite range. The range in which self- 
similarity can be observed (or not observed) changes in 
wave-number space with time, in a way that depends 
both on the initial spectral slope, and on the low-fc and 
high-fc cutoffs in the initial data. 

The paper is organized as follows : in section^ we re- 
call the basic properties of Burgers equation, and give 
a more precise description of the class of initial condi- 
tions we consider. In section IIIII (n > 1 and n < —3) 
we consider the situation when both the velocity and the 
velocity potential are homogeneous Gaussian processes. 
For such initial conditions, we have asymptotically self- 
similar evolution with universal behavior of the spectrum 
E{k) ^ and E{k) ^ k^^ at small and large wavenum- 
ber respectively. For 1 < n < 2, the spectrum, at long, 
but finite time, has also the region |A:|" at very small k 
with time- independent constant, but followed by a re- 
gion fc^ which quickly becomes dominant. In section HVI 
(— 1 < n < 1) we consider the case of homogeneous ve- 
locity potential. In section0(— 3 < n < —1) we consider 
the case of non homogeneous velocity potential. In the 
last two cases the long-time evolution of the spectrum is 
self-similar in some region of (fc, t) plane even we have 
cutoff wavenumber at small and large wavenumber. In 
section lVll we summarize and discuss our results. Details 
of the numerical methods are presented in appendix lVlllI 



II. LARGE-TIME DECAY, SELF-SIMILARITY 
AND BURGERS PHENOMENOLOGY 



We study in this paper the evolution of the velocity 
field, when the initial conditions are random and the ini- 
tial power spectral density is self-similar, that is of the 
form of a power-law Eo{k) ^ \k\" . Let us suppose this is 
the case for a finite interval fci < | fc | < fcu , where fc; and k^ 
are cutoff wavenumbers at large and small scales respec- 
tively - on the infrared and ultraviolet part of the energy 
spectrum. We assume the spectrum to go to zero faster 
than any power-law on either side. We are then inter- 
ested in the plane {k,t), and specifically in the following 
question: where is the behaviour "universal", that is, ex- 



plainable in terms of a few global quantities, and where 
will the specific values of n, fci and fcu play an essential 
role ? 

Fig. n illustrates the results we will show. 

Let us now proceed to explain how Fig.^can be moti- 
vated. From equation Q we can derive an equation for 
the velocity potential ^ and by the Hopf-Cole transfor- 
mation ijj = ivlnO 22, 23], we turn this into a linear 
diffusion equation for an auxiliary field 9. Convolution 
of the initial data for with the standard heat kernel 
gives the the solution of diffusion equation, which in the 
limit v Q may be computed by the method of steepest 
descent. The velocity potential in the limit ^ is then 
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where V'o(a) is the initial potential. The velocity field 
follows by differentiation in x and reads 



v(x,t) = -il)^{x,t) = 



a{x, t) 



t 



(3) 



where a{x,t) is the argument of the maximization in 
for given given x and t. If there are several such a(x, t)'s, 
we are at a shock, where the velocity field is discontuous. 
See e.g. 0,0, IB or op. cit. for an in-depth discussion. 

In this paper we look at initial conditions with energy 
spectrum 



Eoik) = a^\kY'bo(k), 



(4) 



where n is the spectral exponent and 6o(fc) satisfies 
&o(fc) — 1 m some regions of wavenumbers [fci, fcu], while 
going quickly to zero on cither side. The (mean) energy 
is 



E{t) = {v^{x,t)) = E{k,t)dk. (5) 

J — OO 

and the initial energy is denoted 

ul = (vl{x)) = / E„{k)dk. (6) 



It is clear from formula Q above that the solutions 
depend solely from the initial velocity potential. Let us 
introduce the variance of the initial potential (if it ex- 
ists) : 



4 ^ i^oi^)) 



Eojk) 

fc2 



dk. 



(7) 



It is also clear that for a continuous initial velocity field 
the time of first shock formation depends on the initial 
velocity gradients as is, first = —l/Taindxv{x,to). Conse- 
quently, the variance of the initial velocity gradient (if it 
exists) 
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k^Eo{k)dk. (8) 
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FIG. 1: The universal behavior E{k,t) ~ and self-similarity E{k,t) = L^{t)t~'^E{kL{t)) of the energy spectrum 

on the (k — t) plane. The line kL{t) ~ 1/ L{t) is the border between k~^ and |fc|" or k^ asymptotics of the spectrum. The line 
ks{t) is the border between |A;|" and k^ behavior at smal wavenumber. Fig. a - power index n > 2: universal self-similar 
behavior at t S> tni. Fig. Qb power index 1 < n < 2: self-similar behavior only when t 2> fni and k 3> fcs(i). Fig. Qc power 
index —3 < n < 1 : self-similar behavior for t ^ tni. Fig. 0d power index —3 < n < 1 and an infrared cutoff ki : self-similar 
behavior for t-ai and fci ^ ^ fcu. 
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TABLE L Division of domains of spectral exponent n accord- 
ing to second moments of dxvo, vo and i/'o- An entry fci and 
fcu indicates that the corresponding quantity is dominated re- 
spectively by the low-fc and high-fc cutoffs. 



should also be of interest as it determines the typical time 
of first shock formation tg = Since the initial con- 

ditions are scaling only in a finite range, all three charac- 
teristic quantities tTq, cry and exist and are finite. This 
situation we have always in numerical experiments when 
ki is determined by the size of box Lbox and is the 
inverse of the step of discretization. However, depending 
on n, they are dominated by one or the other of the the 
cutoffs. This suggests hence the following first division 
of spectral exponents n (see Table 0. 

From the maximum representation of the solutions to 



n -3 -1 1 

Eit)\ ja^ty' a'k^+'' (Q^)2(n+i)/(3+n) (^^^tyi- 

TABLE II: Division of domains of spectral exponent n ac- 
cording to predicted typical scale and energy of solutions to 
Burgers equation as depending on time t. 

Burgers equation ©, we can introduce the scale L{t), 
proportional to the typical value of \{x — d{x,t)\. For 
large time, balancing the two terms in we have the 
following prediction for the scale L{t) of Burgers turbu- 
lence (see Table inj. Here we take into account that the 
increments of the potential in ^ is {ipo{x) — '0o(O)) ~ 
Q,^(i-n)/2 for n < 1 and is ~ cr^ for n > 1 0|. 
We assume further that ki = for n < 1 and that 
there is some cutoff number k^ for n > 1. In the 
range n < —3 we also need to have fcj > as the so- 
lution of Burgers equation exists only if the potential 
grows slower than quadratically (see the maximum rep- 
resentation Q) and this implies that the spectrum must 
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be shallower than when fc — > 0. From the equa- 
tion 0, we have that at large time between the shocks 
a^shock.m < X < a;shock,m+i the velocity field has an uni- 
versal structure v(x, t) = {x — am)/t and so the energy of 
Burgers turbulence my be estimated as E{t) ^ L?'{t)/t^ 
(see Table For the energy to be finite in the range 
—3 < n < — 1, we require that there is some cutoff 
wavenumber fcj > 0. It has been known for some time 
that the behaviour of L{t) and E{t) in 1 < n has loga- 
rithmic corrections jl^. E^. I25L . 

If indeed Burgers turbulence is characterized by a sin- 
gle scale L{t)^ by dimensional analysis the spectrum takes 
the following self-similar form: 

E{k,t) ^ ^E{kL{t)), (9) 

It is well known that for an initial spectrum with n > 2 
the parametric pumping of energy to the area at small 
fc's leads to the universal quadratic law, E{k,t) ~ fc^, 
and for n < 2 we have conservation of initial spectrum 
E{k, t) = Eo{k) = a^|/c|" at small wavenumber, which is 
the spectral form of principle of "permanence of large ed- 
dies" (PLE) 17, 18] . In Fourier space the self-similarity 
ansatz (O, together with the PLE gives the same rela- 
tions for integral scale and the energy as written above, 
but now in the region —3 < n < 2. Clearly, this argu- 
ment cannot be applied with initial data such that the 
spectral index n > 2, since the later spectrum has now 
a k^ dependence at small k with a time- dependent co- 
efficient. But comparing this with Table 2, where the 
validity of L{t) - (at)2/(3+") is n < 1, we see that the 
region 1 < n < 2 has to be a case apart. In the interval 
1 < n < 2 the self-similarity ansatz is not correct, as was 
shown in [T^ . The reason for this is competition between 
the initial (with constant prefactor a^), and the au- 
tonomously generated fc^ (with prefactor increasing in 
time). If n > 2, the initial spectrum, at low fc, is soon 
overwhelmed by k^ generated by nonlinear interactions 
between harmonics. In this case, hence, the spectrum is 
fully universal, characterized by a single scale L{t), and 
otherwise independent of spectral index n. 

For sufficiently large wave numbers, the spectrum 
should always be dominated by shocks. In one range, 
we should therefore have 

which is equivalent to (O if B{t) ~ L{t)/t^. The am- 
plitude of the small scale part of the spectrum will de- 
crease with time for n > —2 and increase with time for 
-3 < n< -2. 

We note that the spectrum gives only partial informa- 
tion on the solutions to Burgers equation. Indeed, a fc"^ 
tail does not distinguish between discontinuous solutions 
with shocks, and standard Brownian motion, which is 
almost surely continuous. See ^ ^| for other char- 
acteristics of the mass and velocity distribution. 



The rest of the paper will establish the regions in k 
and t for which the above tables is true if we have cutoff 
at large and small wave-numbers. 



III. HOMOGENEOUS VELOCITY AND 
HOMOGENEOUS POTENTIAL; (n > 1) AND 

(n < -3) 



In this chapter we consider the evolution of Burgers 
turbulence in the region n > 1 assuming that both ve- 
locity and potential are homogeneous Gaussian random 
function. It means that we necessarily have some ul- 
traviolet cutoff wavenumber k^^. The function 60 (fc) can 
be characterized by a wavenumber fcu around which lies 
most of the initial energy and which is, in order of mag- 
nitude, the inverse of the initial integral scale Lq. Similar 
situation we have for the spectrum n < —3 and cutoff fci 
at small wave number. 

Most of the results about the energy deca y fo r this re- 
gion have already been obtained by Kida |2J|, but for 
the discrete model of initial condition. He introduced a 
model of discrete independent potential values in adja- 
cent cells, while their relation to the properties of the 
initial conditions (say, the spectrum) were left unspeci- 
fied. For the case of a p.d.f. with a Gaussian tail, he 
obtained the functional form for the correlation function, 
energy spectrum and the log-corrected \/t law for the en- 
ergy decay E{t) ~ t~^a^ ln~^^^ {t/t^\), where however in 
the definition of the nonlinear time t^i was some free pa- 
rameter - the length of cell in the discrete model. In the 
more recent contributions 25], ^ (see also 0) the au- 
thors conjectured the asymptotic existence of a Poisson 
process. This was then proved in |23], showing that, in 
the x-ij} plane, the density of the points is uniform in the 
x-direction and exponential in the ^/'-direction. This per- 
mits the calculation of the one- and two-point p.d.f. 's of 
the velocity |25| and also the full TV-point multiple time 
distributions i\2%. In papers H^,]!^ was also shown 
that the statistical properties of the points of contact 
between the parabola and the initial potential can be 
obtained from the statistical properties of their intersec- 
tions, whose mean number can be calculated using the 
formula of Rice Thus, it is possible to express the 

parameters in the asymptotic formulas in terms of the 
r.m.s. initial potential and velocity. 

In the limit of vanishing viscosity, as the time t tends to 
infinity, the statistical solution becomes self-similar and 
the energy spectrum has the form ||5J). The integral scale 
L{t) and the energy E{t) are given, to leading order, by 

L«).(,,.)V^,n-'"(^), (11) 
£(1) := , (12) 
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where 



tni = L^/cr^, — Lf)/(7y, Lo = cr^/av, 



(13) 



are the nonhnear time and the initial integral scale of 
turbulence. Using this definition we can rewrite in a first 
approximation 

L{t) ^ Loit/Ui)^/^ E{t) ~ Ea{t/U)-\ (14) 

The non-dimensionalized self-similar correlation function 
By{x), which is a function of 5: = x/L{t), is given by 



where for a: > 



1 f 



dz 



(15) 
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.9(2) 



+ g (%^) exp 
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(16) 



(17) 



Our choice of normalization of the energy as E{t) = 
I? (t) /t^ imposes that for the dimensionless spectrum we 
have J E{k)dk = 1. It may be shown that the func- 
tion P{x) is the probability of having no shock within an 
Eulerian interval of length xL{t) Q. 

Note that the properties of the self-similar state are 
universal in so far as they are expressed solely in terms 
of two integral characteristics of the initial spectrum, 
namely the initial r.m.s. potential cr^ and r.m.s. veloc- 
ity (Jv Observe that the spectral exponent does not di- 
rectly enter, in contrast to what happens when n < 1 
(see sections llVIVp . For the dimensionless spectrum 
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By{x) exp{ikx) dx 



ik 



= / xP{x) exp[ikx) dx, (18) 

27r J_oo 



we have the following asymptotic: 

f a I fc2. fc < 1 
E{k) ~ <^ . ' . , k = kL{t) . (19) 
[ a_ k^^, fc > 1 

The k~^ region is the signature of shocks, while the k^ 
region comes due to the parametric pumping of energy 
to the area of small fc. The two constants a-|_ and a_ can 
be computed theoretically as 0] 

poo 

a+ = I/tt / i^P{S:)dx^ 1.078... 
Jo 

a_ = 27r^3/2 ^ 0.359... (20) 



In dimensioned variables, the small-fc region behavior 
of the spectrum is thus 

E{k,t) = a+^r^k^ = A{t)k^, fcL(t) < 1, (21) 



where 



a«,.„,.;ls(±)'i„-v.(_1_ 



(22) 



So, we have a spectrum with an algebraic fc^ region and 
a time- increasing coefficient A{t). 

The situation is more complicated at large but finite 
time 18] . We must now distinguish two cases. When 
rt > 2, the fc^ contribution (|21|1 dominates everywhere 
over the |fc|" contribution and we have a self-similar evo- 
lution in the whole range of wavenumbers (see Fig. ^a)! 
but the "self-similar" time tss from which we have self- 
similar stage of evolution depends on n. In the general 
case, the condition i/ini ^ 1 is not enough for the Pois- 
son approximation to hold and consequently for the ex- 
istence of self-similarity. Let us denote by Acorr a typical 
correlation length for the initial potential, which may be 
greater the initial integral scale Lq- The self-similarity 
occurs when the integral scale of the turbulence L{t) 1)11(1 
is much greater the typical correlation length Acorr ■ This 
leads to the following condition on the "self-similar" time 

tss :i8] 



in 



(23) 



There are instances where (AcorrZ-'^o)^ can be large. Con- 
sider an initial spectrum EQ{k) Q with n » 1 and a 
function 60 (fc) decreasing rather fast when k > k^. In this 
case the initial velocity field is a quasi-monochromatic 
signal with a center wavenumber fcu ~ Lq^ and a width 
Afc ^ [Acorr] ^ fcu. At the early stage of evolution 
tni ^ i ^ ini(Acorr/-^^o) wc havc the saturation of am- 
plitude modulation and the shift of the shocks is much 
smaller then the period of the quasi-monochromatic sig- 
nal. The energy of this signal is approximately the same 
as the energy of the periodic wave: E{t) ~ Lo/12t^ 
[2^. Nevertheless due to the finite width of the ini- 
tial spectrum, we have the generation of a low fre- 
quency component vi{x,t) whose spectrum is well sep- 
arated from the primary harmonic fcu and with the en- 
ergy Ei{t) - £;o(Lo/Acorr)^ < £^0- At t„i(Acorr/£o) < 

t ^ ini(Acorr/io)^ the energy of the low frequency com- 
ponent is larger than the energy of the high-frequency 
quasi-periodic wave, but to the large spatial scale we have 
a relatively small distortion of this component. And only 
at t 3> tss ~ tni(Acorr/-£o)^ do wc havc the self-similar 
regime of evolution. The physical reason for this is a 
strong correlation of the shocks in the early stage of the 
evolution, which prevents the rapid merging of shocks. 
We need to stress that we have a similar situation for a 
spectrum n ^ —3 and a cutoff fci at small wave number. 
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FIG. 2: Evolution of the energy spectrum -with an initial 
spectrum proportional to A;"* at small "wavenumbers k. The 
spatial resolution is A'^ = 2^°. Spectra averaged over 3000 re- 
alizations. The labels correspond to output times t/tni = 
0.033,. . . , up to t/tni = 3.3 X 10^ 




Reduced wavevector, kHt\ 

FIG. 3: Reduced energy spectra E(k,t)i^ j L\^^{t) (n = 4) 
at the same times as a function of reduced -wave number 
kLcxp it) • 

Let us now discuss the results of numerical simulations. 
We use a smooth cutoff of the initial po-wer spectrum 143|l 
"with fcu = iV/16. We consider in all experiments peri- 
odic initial condition, so the infra-red cutoff frequency 
in this case is determined by the size of simulation box 
and fci = 2-K. To check the self-similar ansatz "we consider 
the evolution of energy spectrum E{k,t), of the energy 
E{t) = {v'^{x,t)) and of the integral scale Lexp(t) which 
we can measure from the experimental data as 

Ll^p{t)^{ij'{x,t))/{v'{x,t)) . (24) 

In Fig.Elenergy spectra (averaged over about 3000 re- 
alizations of the random process) are shown at different 
moments of time from t/tn\ — 0.033 to t/t^i — 3.3 x 10®. 
The initial spectrum was at small k. Fig.|21contains re- 
duced energy spectra E{k,t)t^ / Ll^^{t) at the same time 
as a function of reduced wave number kLcxp{t)- We see 
the generation of A{t) k^ with growing amplitude A{t) at 



small k, and k~'^ at large wavenumber. The switch point 
ks{t) between A{t)k'^ and regions of the spectrum 

moves quickly towards the maximum of the spectrum 
and finally the a^k'^ part of the spectrum disappears (see 
Fig.IBi). 

The preservation of the shape of each curve at large 
time t/tn\ > 10 is evident. It is also seen that the to- 
tal energy (the area under the curves) decreases, and so 
does the characteristic wavenumber kL{t) ~ 1/L{t). The 
curves in Fig. |3| have been plotted with respect to the re- 
duced wavenumber k L^xpit) as the quantity Lexp(i) can 
be measured unambiguously from the numerical simula- 
tions. To compare the results with the spectrum the- 
oretical E{k) H18|) . we need to deduce the relation be- 
tween Lcxp{t) and L{t), that we can get from lO 

, r E{k)k'^dk 

Li^jt) = L^t y-Zo ~ ~ ~ 

P J^^E{k)dk 

rOO 

^ L^(t) xP{£)dS: L^(t). (25) 

Jo 

From these, we deduce the experimental values of the 
constants a+ = 1.10 and a_ = 0.37, both slightly larger 
(« 2 %) than the theoretical values. This very small dis- 
crepancy could be due to finite-size effects contaminat- 
ing the measurement of {ijP'{x,t)) at small wavenumbers 
and {v'^{x,t)) at large wavenumbers, and thus of the ex- 
perimental integral scale Lg^pit). 

The asymptotic spectrum is reached rather quickly af- 
ter the non-linear time and the asymptotic formula 
describes the numerical data very well, not only in the 
limits of relatively large and small wave numbers, but also 
at the top, where the spectrum switches between the two 
asymptotes. From Fig. is also seen that the transition 
between the two asyniptotics fc^ and k~^ is rather sharp. 

A. Breakdown of self-similarity; (1 < n < 2) 

The more interesting case is 1 < n < 2, when we have 
breakdown of the self-similarity |0 . The permanence of 
large eddies implies now that, at extremely small k, 

£;(fc,t) ~ a^j/cj", forfc^O. (26) 

This relation holds only in an outer region |fc| <C ks{t) 
where (|26f) dominates over (|21|) . The switching wavenum- 
ber ks{t), obtained by equating (|?T|l and is given by 

Let us define an energy wavenumber k]^(t) — L^^{t), 
which is roughly the wavenumber around which most of 
the kinetic energy resides. From l|ll|l kL{t) ~ (tcr^)~^^^ 
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FIG. 4: Evolution of the energy spectrum with an initial spec- 
trum proportional to |fc|" (n = 1.5) at small wavenumbers k. 
Resolution = 2^". The labels correspond to output times 
ii/ini = 0.18,. . . , up to t/tnl = 1.8 X 10^°. 



(ignoring logarithmic corrections). We then have from 
(|27|l . still ignoring logarithmic corrections: 



kUt) 



(28) 



Hence, the switching wavenumber goes to zero much 
faster than the energy wavenumber, so that the pre- 
served part of the initial spectrum \k\^ becomes rapidly 
irrelevant. Let us also observe that the ratio of the 
energy in the outer region to the total energy, a mea- 
sure of how well the Kida law ((TT|) is satisfied, is equal 
to (t/tni)-3("-^')/("-2) (up to logarithms) and thus be- 
comes very small when t 3> inb unless n is very close 
to unity. Thus, for 1 < n < 2 there is no globally self- 
similar evolution of the energy spectrum at finite time. 
Of course, as n ^ 2 the inner k'^ region overwhelms the 
outer region and as n — > 1 the converse happens, so 
that in both instances global self-similarity tends to be 
reestablished. 

In Fig.0]energy spectra (averaged over ~ 3000 realiza- 
tions of the random process) are shown at different mo- 
ments of time from t/tni = 0.18 to = 1.8 x 10^". For 
the initial spectrum we have n — 1.5. Fig. ^contains re- 
duced energy spectra E{k, t)t'^/L^{t) at the same time as 
a function of reduced wave number kL{t), when we use 
the asymptotic expression (|11() for L{t). We see again 
the generation of A{t) fc^ with growing amplitude A{t) at 
small k and k~^ at large wavenumber. The switch point 
ks{t) between A{t) k^ and a^|fc|" regions of the spectrum 
tends slowly to the origin of the spectrum and finally the 
a^|fc|" part of the spectrum disappears (see Fig. ^b). 

Thus the numerical experiments support the theoret- 
ical prediction that for n > 1 at large time we have the 
self-similar behavior of the spectrum l^. Moreover the 
integral scale of turbulence L(t) ^ll]i and the energy of 
turbulence E(t) are in perfect agreement with theoretical 
predictions, even for 1 < n < 2 where we have breaking 
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FIG. 5: Reduced energy spectra E{k,t)t'^ / (t) at the same 
times as a function of reduced wave number kL{t) using the 
theoretical value L{t) = (i/ini)^/^[log(t/27rini)]"^''''. 




Reduced tiine, // r^,, 

FIG. 6: Evolution of the computed integral scale Lcxp(i) 
for n = 1.5,2,3,4 (dotted lines) compared to the theoreti- 
cal leading-order prediction Lcxp = 1.28 1/(t) Illlll . ll25ll (solid 
line) and the same without the logarithmic correction (dashed 
line) . 



of self-similarity. On Fig. |Hland|7|we plot the evolution of 
Lcxp{t) and E{t) for different n — 1.5,2,3,4 and the the- 
oretical curves, taking into account the relation between 
Lcxp{t) and L{t) (Psjl . showing that the theoretical pre- 
dictions are perfectly reproduced by the simulations. 

In all experiments, we only consider times small 
enough for the characteristic wavenumber k^it) ^ fc;, 
so that we still have many shocks in the box. When this 
condition is broken, it means we have a single shock in 
the simulation domain and we have universal self-similar 
linear decay of the spectrum E{k, t) ^ (kt)^"^. 



IV. HOMOGENEOUS VELOCITY AND 
NON-HOMOGENEOUS POTENTIAL; (-K n< 1) 

We begin with the case —1 < n < 1 when the initial 
potential has homogeneous increments. Many aspects 
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FIG. 7: Evolution of the computed integral energy E{t) for 
n = 1.5,2,3,4 (dotted lines) compared to the theoretical 
leading-order prediction E{t) — L^{t)/t^ 11211 (solid line) and 
the same without the logarithmic correction (dashed line). 

of this case are well understood, thanks in particular to 
Burgers' own work who did consider the case when 
the initial velocity is white noise n = (see also 

The phenomenology is quite simple (see section m . 
Increments A?/'o(i((0) of the initial potential over a dis- 
tance L{t) = X — a{x, t) can be estimated from the square 
root of the structure function S{)^{L) of the potential 129|) 

5o^(a;) = ((^o(x)-Vo(0))'), (29) 

When when n < 1, it grows without bound Sq^{x) ^ 
a^lxj"'^"". For a given position x, the maximum in |2Jl will 
come from those as such that the change in potential is 
comparable to the change in the parabolic term and this 
immediately leads to L{t) - a<^/("+^). (see Table 

When the initial spectrum Q has no cutoff wavenum- 
ber, we can use that scaling and then get immediately 
that the above mentioned expression for the integral scale 
L{t) and energy E{t) (see Table nil are now exact for all 
time starting from zero. For details, see (Section 4) 
and (3^]. But even when the initial spectrum has both 
cutoff waves numbers at large and small scale there is 
some region in (fci) plane where we have self-similar evo- 
lution of the spectrum. 

In the general case for the structure function (|29|l of 
the potential, we have 

So.^{x) = g{x)P^aM^-^ , 

/32 ^ «2 (30) 

r(2 - n) sin ^^i^ 

The properties of the dimensionless function g{x) are de- 
termined by the function h{k) and for h{k) = 1 we have 
g{x) = 1. When the initial spectrum has cutoff wavenum- 
bers ki and ku the function g{x) = 1 in some spatial in- 
terval Lu < a; < Li where Lu ~ 1/fcu and Li ^ l/fci- Let 
us introduce the dimensionless variables 

x = x/L{t), d^a/L{t), d = v/{L{t)/t), (31) 



where L{t) is an integral scale of turbulence. Then the 
"maximum representation" will be rewritten in the form 

v{x,t) = X — a{x,t), (32) 

where a{x,t) is the coordinate at which dimensionless 
function $(i,a, i) achieves its (global) maximum for 
given X and t and 

l>(x,a,t)=^o(a,t)-^^^^. (33) 

here ipQ{d,t) — iljo{aL{t)) {t/L'^{t)^ is the dimensionless 
initial potential with the following structure function 

^ov-(i) =5(ii(0)/3'|ir"" (34) 
Here we define that integral scale L{t) by the relation 

L{t) ^ (ai)2/(3+») (35) 

and so in we use that {atf / L^^+-^'^ (t) = 1. For this 
definition of the integral scale, we have for the dimen- 
sionless spectrum 

{fc" fc <C 1 
, k = kL{t). (36) 
7„fc-2, k > 1 

and for the energy of Burgers turbulence 

E{t) = anL^{t)/t^ = a„a3T^t"TT^ (37) 

where 7„ and an are dimensionless constants, which we 
will determine from the numerical experiments. 

Consider first the case when we have a cutoff wavenum- 
ber fcu only at small scale (fci — 0). For L{t) ^ the 
structure function H35|) my be replaced by 

^n^(i) = /32|5;|i-", (38) 

the function not depending on time and possessing no 
spatial scales in its own. In this case the statistical 
properties of absolute maxima coordinates a do not vary 
with time. The latter phenomenon means that the field 
v{x,t) (|32|l statistical properties determined by time- 
independent statistics of a{x,t) are rendered self-similar 
according 1)32(1 . So at large times when L{t) is also large, 
we have self-similar evolution of the spectrum Eq. @ (see 
Fig. ^c). Alternatively, we could argue that when t is 
so large that the parabolas appearing in ^ have a ra- 
dius of curvature much larger than the typical radius of 
curvature of features in the initial potential, we can plau- 
sibly replace that initial potential by fractional Brownian 
motion of exponent h — (1 — »t.)/2, so that the upper cut- 
off fcu becomes irrelevant. Without loss of generality we 
may assume that this fractional Brownian motion starts 
at the origin for x — Q. This function is then statis- 
tically invariant under the transformation x \x and 
-00 — > A'^^o- It is then elementary, using to prove 
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that a rescaling of the time is (statistically) equivalent 
to a suitable rescaling of x distances and of i^ix, t). This 
implies the above expression for L{t) and E(t) (see Table 

cn). 

For the self-similar initial spectrum with — 1 < ?i < 1 
we have a divergence of the potential in the infrared part 
of the spectrum and divergence of velocity and gradient 
at small scale. So we can not introduce now the initial 
scale on the base of ((T^ . Assuming that we have a cutoff 
wavenumber at small scale, we can in this section use 
some other definition of nonlinear time and initial scale 



(39) 



It is easy to see that this time is equal to the the charac- 
teristic time of shock formation tg. Using this definition 
we can rewrite the expression for L{t) and E{t) in the 
form 



Lit) 

m 



(40) 
(41) 



Assume now that we have also a cutoff wavenumber 
fci at large scale and we have saturation of the potential 
structure function to 2fT^ at a; ^ Lj. In this case in the 
interval when L; ^ L{t) 3> we can replace the dimen- 
sionless structure function (|35|l by (|38|l . It means that in 
some time interval t^^ t ^ ti, where t^, ti are deter- 
mined by the condition L{ti) ~ Li, L{tu) ~ iu, we will 
still have self-similar laws for the integral scale L{t) and 
energy E{t) of turbulence. The energy spectrum E{k,t) 
will have the self-similar behavior on this time in the re- 
gion k > ki (see Fig.^d). For the final state at very large 
times, there are two possible situations. In the first one, 
if we have a finite box of size ibox — Li = 27r/fci and a pe- 
riodic initial perturbation with this period, then at very 
large time we will finally have one triangular wave on the 
period and the energy will be decay as E{t) = Lf/12t^. 
In the non-periodic case with a continuous power spec- 
trum, we will have the generation of a low frequency com- 
ponent E{k,t) ~ A{t) fc^ in region k < ki and finally the 
behavior of the turbulence will be like in the case n > 1 
(see section InH) . 

In numerical simulations we used two models of initial 
spectrum. In the first one we assume that we have the 
self-similar power spectrum in the whole wavenumber 
range from 2tt to Ntt {ki = 2tt, k^ = Ntt), in the second 
we use an infrared cutoff wavenumber 2^ ki ^ k^ — 
Ntt. If the initial spectrum is self-similar inside the in- 
terval ki = 2TT,ku = Ntt then for the nonlinear time t. 
we have 



nl 




(42) 
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FIG. 8: Evolution of the energy spectrum with an initial spec- 
trum proportional to \k\" (n = 0.5) at small wavenumbers k. 
Resolutions N = 2^^ and iV' = 2^° are denoted by filled or 
open symbols respectively. The labels correspond to output 



times ti — 1.5 • 10 



up to t = 1.5 • 10 (in absolute 



time, the reduced time t/t^i depends on the value of the up- 
per cutoff fcu, equal in this case to Ntt or N'tt). 
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FIG. 9: Reduced energy spectra E(k, t) t^l"^ (n = 0.5) at the 
same times as a function of reduced wave number k t'^^'^ . Same 
conditions as in Fig. |H| 



In Fig. |S1 energy spectra E{k,t) are shown at different 

moments of time from t/t^^ — 21 to i/^ni ~ ^-^ ^ ■ 
The initial spectrum was |fc|" with n — 0.5, but two sim- 
ulations were done with different resolutions N = 2^^ 
and N' = 2"^°, so that we have in reality two different 
ultraviolet cutoffs k^ = 2^^tt and k'^ = 2^^tt. Prom the 
figures, one can easily see that in the frequency range 
ki < k < min(fcu,A:^) the spectra are exactly equal to 
each other for all times. Fig. |51 shows reduced energy 
spectra E{k,t)t^/'' as a function of reduced wave num- 
ber k t^l"^ at different times for these two different values 
of N . We see that once again, both spectra display a 
perfectly self-similar evolution and are exactly equal in 
their common range of reduced wavenumbers. Moreover, 
it possible to show that we have not only the conserva- 
tion of the spectrum in presence of high frequency signal 
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FIG. 10: Evolution of the energy spectrum with an initial 
spectrum independent of fc (n = 0) (white noise initial veloc- 

ity). 
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FIG. 11: Reduced energy spectra E{k,t)t''' {n — 0) at the 
same times as a function of reduced wave number kt^^^. 



but also the conservation of the large-scale structures in 
each unique realization [sil l3^ . The measured value of 
the dimensionless constant 71/2 — 1.62 is reported in Ta- 
ble |lll| 

Similarly, in Fig. ^| energy spectra E{k,t) are shown 
at different moments of time from t/t^^ = 1.6 to t/^ni ~ 
5.4 X 10"' for the initial spectrum was a classical white 
noise with n = 0. Fig. 1111 contains reduced energy spec- 
tra E{k, t) t° as a function of reduced wave number k t^^^ 
at five different time. The observed value of the dimen- 
sionless constant 70 = 1.43 is reported in Table ITm 

From the structure of the Burgers equation we see that 
due to nonlinear interaction of harmonics we have always 
the fc^ term at low wave numbers fc, which may be lead- 
ing or sub-leading. The sign of this term depends on the 
initial spectrum. It is evident, that for rt > 2 we have gen- 
eration of new component at small wave number and this 
term increases with time see equation (|22|) . When n < 2, 
this term is subleading, but its growing amplitude can 
make it apparent so that it dominates the dynamics, as 



FIG. 12: Evolution of the difference of the reduced energy 
spectrum with the initial white noise spectrum k" , showing 
the universal k^ term at small ks, which is subleading in this 
case. 




Reduced wavevector, k I ' 

FIG. 13: Evolution of the energy spectrum with an initial 
spectrum independent of fc (n = 0) as a function of reduced 
wave number kt'^^^ with an initial infrared cutoff wavenumber 
at fci = 64 7r (corresponding to 1/32 of the box size). One can 
see that the upper part of the spectrum is unchanged with 
respect to Fig. 1111 

shown in section lTlI Al But when n < 1, this term is com- 
pletely masked by the initial components ~ Ifcl", and the 
only way to show that it is really present is by comput- 
ing the diference of the spectrum with the initial one, 
provided the statistical noise in the simulations is small 
enough. This has been done in Fig. ll2|) . showing this sub- 
leading term ~ fc^, this time with a negative amplitude, 
but displaying also perfect self-similar behavior. 

Another way to display this universal low wavenum- 
ber component, is to introduce some infrared cut- 
off in the initial condition. In the simulation shown in 
Fig. ^1 we consider the case of white noise initial spec- 
trum n ~ but with an infrared cutoff at fcj = 64 tt. 
For the first two displayed times, we have a self-similar 
evolution of the spectrum in the wave-number range 
fci < fc < fcu, but with the generation of a component 
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FIG. 14: Evolution of the computed integral energy Eit) for 
n = 1,0.5,0,-0.5. 



n 


-2.5 


-2.0 


-1.5 


-1.0 


-0.5 


0.0 


0.5 


1.0 






0.94 


0.97 


1.09 


1.26 


1.43 


1.62 


1.93 












0.43 


0.32 


0.32 


0.37 



TABLE III: Measured values of the universal constants 7„ 
and a„ characterizing the fe~^ spectrum shock tail and law 
of energy decay respectively. Values marked — don't exist or 
couldn't be measured. 

spectrum E{k,t) ^ A{t) in the region k < ki. At the 
time = 1.6-10'^, L(t) ~ Lj reaches the lower cutoff fcj, 
and the spectrum then becomes equivalent and evolve in 
time as in section lllll where we have self-similar evolution 
with universal behavior of the spectrum E{k) ~ fc^ and 
E{k) ~ k~'^ at small and large wavenumber respectively. 
It is interesting to note that all parts of the spectrum fc^ , 
k^ and fc~^ evolve in a self-similar way for intermediate 
times ijji < t < ti. 

On Fig.^Jthe evolution of the energy E(t) for different 
n — 1, 0.5, 0, —0.5 is plotted. The law of decay is in 
good agreement with the theoretical prediction H41|l for 
all time. The constants a„ have been measured and are 
shown in Table ITTll 



V. NON-HOMOGENEOUS VELOCITY AND 
NON-HOMOGENEOUS POTENTIAL 
POTENTIAL; (-3 < n < -1) 

For the self-similar initial spectrum with — 3 < n < — 1 
we have a divergence of the potential and the velocity 
in the infrared part of the spectrum, and of the gradient 
in the ultraviolet part. Assuming that we have a cut- 
off wavenumber at small scale ki, we can still use the 
definition of the nonlinear time through the gradient of 
velocity t^^ = l/aq which is equal to the the charac- 
teristic time of shock formation ts- Due to the divergence 
of the energy in the infrared part of the spectrum, the 
dissipation in the shocks doesn't lead to a finite value 
of the energy at any time, if there is no cutoff fc;. Nev- 




10° lo' 10* 10° 



Waveveclor, k 

FIG. 15: Evolution of the energy spectrum with an initial 
spectrum proportional to k~^. One can see that the spectrum 
doesn't change at all until the non-linear time of the smallest 
wavenumber component is reached (near t « 0.15). 



ertheless we can still introduce the integral scale of tur- 
bulence L{t) showing the region when initial power law 
spectrum E(k, t) ^ |fc|" transforms to the universal spec- 
trum E{k,t) ^ . For the spectral form of principle 
of "permanence of large eddies" , we still have that inte- 
gral scale grows according ((23 • From equation ((Tn|l . we 
see that the amplitude of small scale part of the spec- 
trum decreases for n > —2 and increases with time for 
-3 < n < -2. 

Let as start for the special case of initial spectrum 
with the critical index n — —2 when L{t) ^ (at)'^. From 
equation we see that the spectrum does not change 
in time. In Fig. 1151 energy spectra E{k,t) are shown at 
different moments of time from t/t^^ = 0.005 to = 
5. And we really see that it is only when ^ 1 

that the spectrum begins to change, simply decaying in 
amplitude without changing its shape. But even if it is 
not apparent in the spectrum, we have an evolution in 
time of each realization and of other statistical properties, 
like shock probability distribution or higher moments of 
the velocity. 

On Fig. ^1 the evolution of a realization of a velocity 
field with a k""^ spectrum is plotted at different times. It 
is easy to see, that even if the spectrum does not change, 
the characteristic distance between the shocks which is 
proportional to the scale L{t) increases with time. The 
signal switches continuously from a Gaussian Brownian 
motion with random phases to a triangular wave with 
aligned phases, but these signals have the same power 
spectrum, and the amplitude of the spectrum is pre- 
served by the Burgers evolution. The evolution of the 
spectrum in numerical simulations starts and the energy 
begins to decay at the time tni,i when the integral scale 
of turbulence L{t) reaches the size of the box Lbox- This 
can also be seen in Fig. El showing the spectra rescaled 
with time for Gaussian Brownian motion initial condi- 
tions. One can see that the initial fc~^ spectrum visi- 
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FIG. 16: Evolution of the velocity field for one particular 
realization with an initial spectrum k~^. Upper plot is the 
Gaussian Brownian motion initial condition att — and lower 
plot is the velocity field at t = 0.15, just before the decay of 
the smallest wavenumber component. Both signals have the 
same power spectrum. Note that the largest scales are non- 
hnearly distorted, but their amplitude is globally preserved, 
illustrating the principle of persistence of large eddies. 
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FIG. 17: Reduced energy spectra E{k, t) as a function 
of the reduced wavenumber k L{t) = kt^ for an initial spec- 
trum Note that the spectrum rescaling works even if this 
case with a factor decreasing with time, and allows the sepa- 
ration of the k~^ spectra in difi'erents zones coming from the 
initial conditions (k L(t) < 1) of from the shocks (kL{t) > 1). 



ble for k L{t) < 1 is continuated "with the spectrum 
of the shocks at late times for kL{t) > 1, but that the 
rescaling allows the separation of both. For a given finite- 
size realization of the signal "with both lo"wer and upper 
wavenumber cutoffs fc; and fcu, the spectrum will slide 
along the curve in Fig. 1171 until kiL{t) > 1, that is the 
rescaled wavenumbers k L{t) > 1 are larger than 1 for 
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FIG. 18: Evolution of the energy spectrum with an initial 
spectrum proportional to One can see that the am- 

plitude of the high-frequency components increase with time, 
the spectrum transforming to a k~^ tail, whose en- 

ergy really comes from the smallest wavenumber components. 
These have a very slightly decreasing amplitude, so that the 
global energy is nearly constant, but slightly decays in time. 



all fcs, at which time the spectrum will begin to decay 
linearly (as shown in Fig. E|for t/t^y > 27). 

In fig. ^1 the energy spectra E{k, t) with an initial 
power law spectrum E(k,t) ~ |A;|" with n ~ —2.5 
is shown at different moments of time from t/t^^ — 
6.9 X 10^^ to i/ij^j = 6.9. We see the generation of 
universal k^^ tail which amplitude increases with time 
and switching point between k~'^'^ and parts of the 
spectrum moved to the small wavenumbers. The value of 
dimensionless constant 7„ in the dimensionless spectrum 
in (|36|l are shown in Table IIIII Note that the constant 
for n — —2.5 has not been measured as the energy begins 
to decay immediately as soon as the tail is estab- 
lished. This effect also perturbs the measurement of 7„ 
for n < — 1, because in a finite-size system, some (small) 
dissipation occurs as soon as shocks are present, and so 
the amplitude of the k~'^ tail is diminished. This explains 
why we measure for instance an amplitude 7_2 = 0.94 at 
large times, even if 7_2 = 1.0 should be observed. 

We can also introduce the other nonlinear time is a 
time of nonlinear decay of first harmonic tn\,i ~ ki/Ai. 
For — 3 < n < 1, we have no significant decay of the 
energy of turbulence until t > tni,i. The constancy of 
the energy until much later than t^^ is also evidenced in 
Fig. El for n = —1.5,-2,-2.5. As the energy doesn't 
decay until ini.i, the "constants" a„ don't exist for n < 
— 1 and are thus not shown in Table ITTll 



VI. DISCUSSION 

In this work we have reconsidered again the classical 
problem of the spectral properties of solutions of Burgers' 
equation for long times, when the initial velocity and 
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FIG. 19: Evolution of the computed integral energy E{t) for 
n = —1.5, —2, 0, —2.5 as a function of t/t^i. 

velocity potential are stationary Gaussian processes. We 
have shown in greater detail how the self-similar (and not 
self-similar) regimes are realized with initial conditions 
that are only self-similar over a finite range. The range 
in which self-similarity can be observed (or not observed) 
changes in wave-number space with time, in a way that 
depends both on the initial spectral slope, and on the 
low-fc and high-A: cutoffs in the initial data. 

Depending on the statistical properties of initial veloc- 
ity and potential on can introduce the following regions 
on n-axis. They are: homogeneous velocity and homoge- 
neous potential: {n > 1) with subinterval (1 < n < 2); 
homogeneous velocity and and non-homogeneous poten- 
tial: (—1 < n < 1); non-homogeneous velocity and non- 
homogeneous potential: (—3 < n < —1) with some criti- 
cal point n = —2. The common properties of turbulence 
is the self-similar behavior, determined by only one scale 
L{t) - the integral scale, for a range of time and wave 
numbers even in the presence of high- or low-frequency 
cutoffs. But the type of self-similarity is different in 
different region on "n-axis" and is determined by the 
properties of initial potential. High-resolution numerical 
simulations have been performed confirming both scaling 
predictions and analytical asymptotic theory. 

For the case (n > 1) we have verified by numerical ex- 
periments the asymptotic theor y d erived previously by 
several groups [ll H IH IH II^ . The main results 
are : at very large times the spectrum tends to a limiting 
shape, proportional to fc^ at small wave numbers, and 
to at large wave numbers, such that evolution shape 
is determined by peak wave number, fci,(i) ~ l/L{t). 
Due to the merging of the shocks the integral scale in- 
creases with time L{t) (ta^)^^"^ In""^^** (tfT^/27rcr^). So 
asymptotically the evolution of Burgers turbulence and, 
in particular, the law of energy decay is determined only 
by the variance of potential tr^. 

For large, but finite time, we have breakdown of self- 
similarity of the spectrum when 1 < n < 2. The spec- 
trum then has three scaling regions: first, a |/c|" region 



at very small k < ks{t) with a time-independent con- 
stant, second a k^ region at intermediate wavenumbers 
ks{t) < k < kL{t) with increasing amplitude, and, fi- 
nally, the usual k^^ region at fc > k^lt). The relative 
part of the spectrum with |fc|" region decreases with time 

h{t)/kL{t) ^ (t/tnl)"'""''/'*'""^- 

In the case of finite viscosity, if one introduces an 
instantaneous Reynolds number Re(t) ~ L(t)E^/^{t)/i' 
based on viscosity, the typical velocity and typical spatial 
scale at time t, it means that Re(i) ~ Keo{logt/tni)~^/^ . 
Within dimensional estimates, the Reynolds number 
would be constant in time. On a practical level, we have 
thus established that the inviscid approximation is not 
valid for arbitrary long times. After a time, which is 
very long if the initial Reynolds number is large, about 
iniexp(ReQ), the viscous term in becomes comparable 
to the inertial term everywhere, and can not be neglected. 

When n is less than one, the large-scale part of the 
spectrum E(k) = a^|fc|" is preserved in time and the 
global evolution is self-similar, so that scaling arguments 
perfectly predict the behavior in time of the integral scale 
L{t) = (ai)2/(3+»). For -1 < n < 1, the energy also 
decay as a power law E{t) - Q,4/(3+n)i-2(«+i)/(3+n)_ 
case of finite viscosity the increasing of the integral scale 
is faster the decay of the energy and we have for the 
Reynolds number Re{t) - i(i-")/(3+n)^ i.e. the Reynolds 
number increases with time and the shape of the wave 
becomes more and more nonlinear. This last point is true 
only in the case when we have not cutoff wavenumber 
at large scale. In numerical simulations in a finite box, 
the final behavior will always be the linearly decaying 
sinusoidal wave with period equal to the size of the box. 
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VIII. APPENDIX. NUMERICAL WORK 



A. Normalizations of initial spectrum 

We use the following smooth cutoff of the initial power 
spectrum : 

_ fc^ 

Eo{k) = a^\k\"e ^. (43) 
The variances of the velocity, the velocity gradient and 
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potential can be computed to be: 
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1 




2 




n 
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3 




2 




n 




1 


2 



l.n+1 
'^0 ' 



l.n+3 
'^0 ' 



n-1 



(44) 
(45) 
(46) 



where T{x) is the gamma- function. 

From this equation it is easy to see the critical points of 
n connected with the divergence at small wave numbers. 



B. Generation of initial conditions 

The scale of the box in the numerical experiment is 
taken as unit of space, so that wavenumbers k go from 
27T to Ntt, where TV is the number of points used in the 
simulation, typically ranging from N = 2^^ ~ 3 x 10^ to 
N = 2^" ~ 10^ in our simulations. The amplitude of the 
spectrum (143(1 was simply taken as a = 1. Fourier com- 
ponents of a Gaussian process are independent Gaussian 
variables. We therefore synthesize the initial potential of 
the velocity by first generating random Fourier compo- 
nents Ofc distributed according to 



p(afe) 



: exp 



where 



a1 — E^{k)dk. 



And we use the well known relation between the power 
spectra of the process and of its integral 



E^{k) = fc-2^o(fc) 



(47) 



where the form of _Eo(fc) is chosen with a smooth cutoff 
at large k according to H43|l . 

By inverse Fourier transforming the components, we 
obtain the initial potential in real space, from which we 
can obtain the potential at any time using the Legendre 
Transform Q (see 9, 33] )■ Repeating the whole process 
many times with different realizations of ipo{x), we sam- 
ple the desired ensemble of Gaussian initial conditions. 



C. Fast Legendre Transforms 



In numerical simulations the initial data are always 
generated as a discrete set of N points. It could be as- 
sumed naively that the number of operations necessary 
to compute the maximization Q for all values of x scales 
as 0{N'^). It may however be shown, using (0) that 
a{x) is a nondecreasing function of x. The number of 
operations needed in an ordered search therefore scales 
as 0{N \og2 N) when using the so-called Fast Legendre 
Transform procedure [3, U3, ■ 
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